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ABSTRACT 


The subject of this thesis is an investigation of the effect of using 
the lambda-sigma jump process in the acoustic detection component of 
APAIR. A computer simulation was developed which is similar to the 
sonobuoy field versus submarine engagement model found in APAIR, the 
Navy’s general ASW model. This simulation was then modified to 
incorporate the lambda-sigma jump process and the effect of this 
modification is discussed. In order to check the structural validity of 
the simulation models, results that were obtained by using them are 
compared to results that were obtained by using an analytical model called 


the random search model. 
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I. INTRODUCTION 


The number of sonobuoys needed to detect a submarine with a given 
probability can be estimated using a computer simulation. Various passive 
acoustic detection models that have been developed over the years can be 
used to form the basis of this simulation. One of these detection models 
is the stochastic signal excess model that is part of the sonobuoy field 
versus submarine section of APAIR, the Navy’s general air ASW model. 
Another detection model that has been extensively used is one that 
incorporates a lambda-sigma Jump process. 

The purpose of this thesis is to investigate the effect of using the 
lambda-sigma jump process in the passive acoustic detection component of 
APAIR. To do this, a computer simulation, called Model A, was developed 
that incorporates the detection model found in the APAIR model. Model A 
was then modified by incorporating the lambda-sigma jump process in order 
to obtain a second computer simulation, called Model B. This modification 
was done in order to determine if that change would affect the estimate 
of the number of sonobuoys required to gain a detection with a given 
probability. 

Chapter 2 provides a description of the two simulation models. Data 
generated by the two simulation models are compared in Chapter 3 with data 
obtained using an analytical model called the random search model. This 
was done in order to provide a check on the simulation’s structural 


vauidity 


Il. TRE He@BiS 


Three models are discussed in this chapter. The first two are 
Simulation models that are based on a stochastic signal excess acoustic 
detection model. The first model is called Model A and the second model 
is called Model B. The third is an analytical model called the random 


search model. 


A. MODEL A 

Model A is based on the following assumptions: At the start of a 
search, all the sonobuoys are in the water and operating and the submarine 
iS in a submarine probability area, or SPA. The size of the SPA is 
4900 nm° and it is assumed to be Square in shape. The submarine’s initial 
location is uniformly distributed over the SPA, its speed is constant at 
5 kts, its course is constant and is equally likely to have any value 
between 0 and 360 degrees and its depth is 400 feet. The sonobuoy’s 
receiver depth is 400 ft. 

The environmental data for the model, such as propagation loss and 
ambient noise, are based on an area in the northwest part of the pacific 
ocean. Since some of the sonobuoys that are currently in use or are being 
developed may not have the ability to achieve convergence zone detections, 
it is assumed that detection can only be made out to a maximum of 10 nm 


in order to eliminate the possibility of this type of detection. 


aid eeeesseiceune diiierence Detween the signal-to-noise ratio in 
decibels and the detection threshold. If the signal excess is greater 
than or equal to zero, then the probability of detection is greater than 
or equal to a minimum acceptable value for a single observation. In this 
thesis, the minimum acceptable value is .5 and recognition differential 
is the term that is used to refer to detection threshold. For the model 
and for any time t, signal excess is a random variable that is defined as 
follows: 

ee Acct) = SE(t) + X(t). 
In equation 1, X.-(t) is the signal excess, SE(t) is the mean of the signal 
excess and X(t) is a random variable that determines the stochastic 
Bimeac eer OT the Signal excess. Since SE(t) is the mean of X,-(t), the 
mean of X(t) is equal to zero and the standard deviation of X(t) is equa] 
to the standard deviation of X.-(t). The signal excess can also be defined 
as follows: 

By re Xe (8) Ut) - (xy (t) - KEY) - p(t). 
The subscripts have the following meaning: SE indicates signal excess, SL 
indicates the target’s source level, TL indicates transmission loss, NL 
indicates noise level, DI indicates directivity index and RD indicates 
recognition differential. It is usually assumed that al] the random 
variables on the right hand side of Equation 2 are normally distributed 
and statistically independent and thus X.-(t) is normally distributed and 
has a variance equal to the sum of the variances of the random variables 


on the right hand side. 


The expected value of X.-(t) as defined by (2) is: 

(3) SE = SL = lL - (ies 0! | ae 
which can be interpreted as the passive sonar equation. The passive sonar 
equation represented by (3) can also be written in terms of the figure of 
merit, or FOM, as: 

(4) SE = FOM - TL 
where 

(5) FOM = SL - (NL - DI) - RD 
(see appendix A for a discussion of the passive sonar equation). The 
detection model that is the basis for the computer simulation is similar 
to the one used in a model developed by Vitro called RADS, an acronym for 
Rapid Acoustic Detection Simulation [Ref. 1]. The RADS detection model 
in turn, is similar to the sonobuoy field versus submarine acoustic 
detection model used in the APAIR model, the Navy’s general air ASW 
engagement model. RADS was developed to specifically simulate a sonobuoy 
field’s detection capability. Faster results can be obtained by using 
RADS as opposed to APAIR when one is interested only in sonobuoy field 
detection capability. Unlike the RADS model, Model A does not take into 
account aircraft motion and only simulates the target and the sonobuoy 
field. Also, it is not capable of simulating a transiting battle group 
performing barrier operations against a submarine. However, like RADS it 
is a monte carlo simulation of a sonobuoy field containing a single 
Submarine. 

In Model A, a submarine is treated as a point sound source which is 


not aspect dependent. Up to three separate source frequencies may be 


specified at three separate source levels. The source levels which are 
determined off-line by the user, depend on the submarine type, speed, 
depth, and the operating conditions. For this thesis, only one frequency 
was used. 

The simulation run time or game time is user controlled. For this 
thesis, it was 1200 minutes. This value was chosen since, for a submarine 
travelling at a constant speed of 5 knots, the maximum time the submarine 
can remain within a 4900 nm Square sonobuoy field is approximately 1200 
minutes. The sonobuoys are placed in the search area such that there is 
equal spacing between each of the buoys and between the outer buoys and 
the edge of the SPA. It is assumed that all the sonobuoys operate for the 
entire 1200 minute run time. 

The sonobuoy and acoustic processor characteristics are controlled by 
the user by inputting RD, and DI. 

The primary features of the ocean acoustic environment which influence 
passive sonobuoys are ambient noise and transmission loss. These 
parameters are supplied by the user as part of the input data set. Three 
values of ambient noise can be input, one for each source frequency. 
Transmission loss is input in the form of tables, one table for each 
source frequency. The table consists of transmission loss in dB for every 
one-half nautical mile increment out to a user controlled maximum limit. 

The detection process is simulated by computing the signal excess for 
each source frequency at each sonobuoy during every time step in the 
Simulation. The time step increment is controlled by the user and remains 


constant. At each time step of the simulation, the signal excess for each 


source frequency is determined as follows: First, the range between each 
Sonobuoy and the submarine is computed. Based on this range, an 
interpolated value of transmission loss is calculated. Then, the signal 
excess is determined by using (1) where: 

(6) SEG) eaSiecaNib - TUGt ReDim 
and: 

(Ze) Xi st) aXe aia ety Rie 

This model simulates signal excess fluctuation by separating the total 
fluctuation into three individual and independent components: (1) long 
term or day to day fluctuation; (2) short term or minute to minute 
fluctuation; and (3) buoy to buoy fluctuation. 

Long term fluctuation is caused by variation between operators, 
variation between average source levels of a particular submarine type and 
Variation in ocean acoustic conditions caused by seasonal changes and day 
to day changes [Ref. 2]. X, represents the long term fluctuation and is 
a random variable whose values are determined by drawing from a normal 
distribution with a mean of zero and a user supplied standard deviation. 
Once the value of X, is obtained at the beginning of the simulation it 
remains constant for a full replication and the same value is used for 
every signal excess calculation at every sonobuoy. This implies that the 
model assumes complete dependence of the sonobuoys with respect to long 
term fluctuation. 

Short term fluctuation is caused by changes in target aspect, changes 
in operator alertness, and variations in ocean acoustic conditions caused 


by shipping traffic and wind [Ref. 2]. X,(t) represents the short term 


fluctuation. It is a random variable whose value is determined several 
times during a replication. The standard deviation and sampling interval 
or time between random draws of X,(t) are controlled by the user. At each 
sampling interval, a random draw is made from a normal distribution with 
a mean of zero and a user supplied standard deviation. The sampling 
interval of random draws of X,(t) was set at 30 minutes for this thesis 
[Ref. 3]. The values used in the signal excess equation for X,(t) are 
interpolated values between the values determined by successive draws of 
X,(t). Just as with X,, the same value for X,(t) is used to compute signal 
excess at every sonobuoy. Again, this implies complete dependence between 
sonobuoys for short term fluctuation. 

The final term X, represents buoy to buoy fluctuation. This type of 
fluctuation could be caused by variations in the manufacturing process of 
the sonobuoys. At the beginning of each replication, a random draw is 
made for each sonobuoy from a normal distribution with a mean of zero and 
a user supplied standard deviation. This X, value is determined separately 
for each sonobuoy and is held constant for an entire replication. Once 
the three separate fluctuation terms are computed they are combined to 
form X(t) and added to the sonar equation to obtain a signal excess value. 

As stated earlier, the user must input values for the standard 


deviation of the random variables X,, X,(t), and X,. These values are not 


readily available. 


A sensitivity analysis was conducted on these standard deviation terms, 
to determine if varying their values would have a significant effect on 
the probability of detection by using the following assumptions: (1) the 
Standard deviation of signal excess is 8 dB [Ref. 4]; and (2) the 


individual standard deviations can be varied subject to: 


(8) oman Or oo, 

Where o,, is the standard deviation of signal excess, o, is the standard 
deviation of the long term fluctuation, o, is the standard deviation of the 
short term fluctuation and o, is the standard deviation of the buoy to buoy 
fluctuation. 

The signal excess based on a standard integration time is computed at 
each time step. Detection occurs on the time step if a moving time 
average of the signal excess is above one of a set of user controlled 
Signal excess decrements. These decrements account for the increase in 
RD at the initial time steps that result from integration times less than 
the standard integration time. The moving average is (S, + Si). +... 


+ Sacieay is where j = 1, 2, 3,... is the time step index and k ee 


tl 


3,...,;J, if J < M and k = M if j >= M where M is the number of samples 
in the standard integration time. The signal excess decrements are input 


in the form of a table of values DEC(i). 


For every sonobuoy and frequency, detection will occur when one of the 
following expressions is satisfied: 

S, >= DEC(1) 

Sees) © >= DEC(2) 


ee et Pg) >= DEC(M). 


An example of a DEC table for M = 5 is as follows: 


DEC(1) = 3.5 
DEC(2) = 2.5 
DEC(3) = 1.5 
DEC(4) = 0.5 
DEC(5) = 0.0. 


This process of time integration of an acoustic signal is used in a model 
described by Forrest [Ref. 2]. Also, a continuous time model is described 
by McCabe and Belkin [Ref. 5]. It can be seen from this table that if the 
initial signal excess sample is 4 dB a detection will occur at the initial 
time step. If the initial Sample is 3.2 dB and the second sample is 2.2 
dB, neither sample is sufficient to make a detection. However, the time 
average of 2./ dB is sufficient to make a detection based on the DEC(2) 
value of 2.5 dB. This simulates the decrease in the detection threshold 
with increasing integration time. The DEC values were obtained using a 


detection model described by Forrest [Ref. 2:pp. 27-29]. 


In this model, the probability of false alarm and the probability of 
detection can be defined as: 

yep, = hee 
and: 

(10) py = &(-v + d!/”?) 
where p, represents the probability of false alarm, py represents the 
probability of detection, $¢ symbolizes the standard normal cumulative 
distribution function, vis a threshold value and d is the detection 
index. The detection index is defined as: 

(11) d = t(BW)(S/N)° 
where t represents time, BW represents bandwidth, and S/N represents the 
Signal-to-noise ratio. Since detection threshold is defined as: 

(12) DT = 10 log(S/N) 
by using equation (11), equation (12) can be rewritten as: 

(13). Di se5alog (due 
Assuming a p- of 0.001 and a p, of 0.500 and using equations (9) and (10) 
the detection index, d, is computed to be equal to 9°55]. Since pam 
0.500, RD can be used vice DT. Assuming a BW of 1] Hz, a solution for RD 
can be obtained for various values of t. For this thesis, the integration 
time and the signal excess sampling interval were 12 minutes and two 
minutes respectively. Table 2.1 shows the values of RD for the various 
values of t. These 6 RD values were then used for the signal excess 


decrement values in the simulation. 


Table: 2.1 
RD values for corresponding 
values of t 


8, Gels: 
1.890 
1.010 
Wp SiS38, 
S000 
-0.496 
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Either a one, two, or three-source frequency (line) detection may be 
specified by the user. If, for example, a two-line detection is required, 
then one sonobuoy must be in contact on two frequencies in order for a 
detection to occur. One frequency held by one sonobuoy and another held 
by a second sonobuoy does not constitute a two-line detection. As noted 
above, only one frequency was used and therefore a one-line detection was 


specified. 


B. MODEL B 

A second simulation model called Model B was obtained by modifying the 
way Model A calculates X(t). Model B uses a lambda-sigma jump process to 
calculate X(t). For a description of the lambda-sigma jump model see 
Reference 5. This modification was made in order to determine if the 
method by which the two simulation models calculate X(t) would have a 


Significant effect on the models estimate of the number of sonobuoys 


1] 


required to gain a detection with a given probability. The standard 
deviation of X(t), for Model B, was assumed to be 8 dB and the mean time 


between jumps was 30 minutes. 


C. THE ANALYTICAL MODEL 

This model which is a form of the random search model is based on a 
deterministic detection model. The model has been used in making 
expected-value calculations of the probability of detection as a function 
of time for a submarine moving within a sonobuoy field [Ref. 6:pp. 28-29]. 
It was used in the investigation in order to provide a check on the 
Simulation’s structural accuracy. 

The performance of each sonobuoy is characterized by an effective 
detection range which depends on the submarines radiated signal level, the 
environmental parameters and the searcher’s signal processing capability. 
The sonobuoys are assumed to be systematically distributed throughout the 
search area such that their coverage is non-overlapping. 

The following parameters appear in the model: 

Q: effective sonobuoy detection radius. 

N: number of sonobuoys in the field. 

V: submarine speed. 

S: average sonobuoy spacing. 

T: search time. 

It 1s assumed that Q is equal to the median detection range or MDR 


which was obtained from the transmission loss data used in the simulation. 


yz 


For this thesis T = 6.56 hours and is the average time the submarine 
remains within the sonobuoy field. 

The submarine can be detected in two ways: First, it can be detected 
at the beginning of the search if its initial position is such that it is 
within the detection range of one of the sonobuoys. Second, it can be 
detected if the submarine’s initial position is such that it is outside 
the detection range of all the sonobuoys but subsequently it moves into 
detection range of one of the sonobuoys. This leads to the following 


expression for the probability of detection [Ref. 6:pp. 28-29]: 


P= 1 - ((1 - (02/52) Je ~((2QVTI/(S°CL = (Q°/S°))))). 
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lil. RESUETS 


This chapter presents the results obtained from Model A, Model B and 
the analytical model. Each of the simulations was done using FOM values 
of 55 dB, 60 dB, 65 dB and 70 dB. The results of the simulations were 
then compared with the results obtained using the analytical model. In 
Model A, a sensitivity analysis was done of the three standard deviation 
terms of the random variables X,, X,(t), and X,, to determine their effect 
ONaua 

Probability of detection for an encounter between the submarine and 


the sonobuoy field was estimated using the statistic: 


’ number of detections 
replications 


One hundred replications were used. Ideally, a larger number of 
replications would have been used in order to obtain a higher degree of 
Statistical validity. However, one hundred replications is a compromise 
that was made in order to limit computer processing time. Pseudorandom 
numbers were generated using the LLRANDOMII package as installed on the 
IBM 3033 at the Naval Postgraduate School [Ref. 7]. | 
Figures 3.1 through 3.4 provide a graphical comparison of the 
probability of detection as a function of the number of sonobuoys for each 
of the models. The data corresponding to these curves is summarized in 


Tables 3.1 to 3.4. The values used for o,, 06., and o, in Model A were 4.6 
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dB, 4.6 dB and 4.6 dB respectively, corresponding to a standard deviation 
of 8 dB for X(t). For Model B the standard deviation of X(t) was 8 dB. 

A sensitivity analysis was performed on the three standard deviation 
terms and the results are graphically presented in Figure 3.5 with the 
corresponding data in table 3.5. Figure 3.5 was generated by running 
Model A, using four combinations of the three standard deviation terms at 
an FOM of 60 dB. For three of the combinations, each standard deviation 
term, 6,, 0, and o, designated as L, S and B in the graph, was set at 7.0 
dB while the remaining two were set at 2.7 dB. The fourth case had the 
three standard deviation terms equal to 4.6 dB. 

Smoothed curves were superimposed over the data generated by the two 
Simulation models, in order to gain a better impression of the pattern of 
dependence of P, on the number of sonobuoys. A procedure known as locally 
weighted regression or LOWESS was used to generate the smoothed curves 
[Ref. 8:pp. 94-104]. 

An inspection of Figures 3.1] through 3.4 indicates that there were 
differences between the three models. The relationships between the three 
curves seem to be independent of FOM. As FOM decreases the number of 
Sonobuoys required to obtain a given probability of detection increases, 
for all three models. The curve generated by the analytical model lies 
between the two curves generated by Model A and Model B and the curve 
generated by Model B lies well above the curve generated by Model A, for 
al] values of FOM. 

The difference between Model A and Model B is in the way they model 


the stochastic component of signal excess and this difference affected the 
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estimate of the required number of sonobuoys to gain a detection with a 
given probability. The curve generated by Model B which corresponds to 
a lambda-sigma jump signal excess fluctuation process lies well above the 
curve generated by Model A. So, for example in estimating the number of 
sonobuoys required to obtain a 50 % probability of detection of the 
Submarine at an FOM of 55 dB, Model B estimates it would take 
approximately 75 sonobuoys, Model A 150 sonobuoys and the analytical model 
125 sonobuoys. 

An inspection of Figure 3.5 indicates that changing the values of the 
Standard deviation terms causes an effect on the simulation models 
estimation of the number of buoys required to gain a detection with a 
given probability. The curves start out fairly close together and then 
Separate as the number of sonobuoys increase. The curve generated by the 
case when the long term fluctuation term was dominant lies below the 
curves generated by the other three cases as the number of sonobuoys 
increase past 100. 

To determine if changing the values of the standard deviation terms 
had a statistically significant effect on the estimate of the number of 
sonobuoys required to gain a detection, a contingency table was used and 
Appendix D presents a discussion of this method. The conclusion using the 
contingency table was that there was a dependence between the probability 


of detection and the three standard deviation terms, 0,, oO, and o,. 
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Finally, a few comments should be made concerning the results. For 
this scenario it seems the analytical model might be as capable and 
certainly less costly with respect to the two simulation models in 
estimating the number of sonobuoys required to gain a detection. 
Model B consistently produced higher values of P, than Model A. This could 
be due to the fact that with Model B there is greater independence between 
Successive values of X(t) whereas with Model A that is not the case. As 
the independence between successive values of X(t) increases so should the 


probability of detection. 
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Table 3.1: Probability of Detection vs. Number of Sonobuoys 
for a Figure of Merit of 55 


Number of Probability of Detection 
Sonobuoys 
Analytical Model B 
4 0.02 O,. 0.06 
g 0.05 0. ORO 
16 0.08 Oz Oars 
Z5 013 OF 0.24 
36 0.18 Oz 0.18 
49 0.24 OF 0.46 
64 0.30 0. 0.54 
8] 0.36 OF 0.49 
100 0.42 0. 0.63 
WZ) O50 OF Oz 2 
144 Oe 5 OF 0.81 
169 0.61 0. 0.68 
196 0.66 OR 0.96 
Va As) Oey On 0.9] 
256 6.76 OF 0.98 
289 0.8] aly 0.95 
324 0.84 0. 1.00 
36] Ore BF On 1.00 
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Figure 3.2: 
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Table 3.2: Probability of Detection vs. Number of Sonobuoys 
for a Figure of Merit of 60 


Number of Probability of Detection 
Sonobuoys 

Analytical Model B 

4 0.04 0. 0.08 

9 0.08 0. Oe 
16 0.14 OF 0325 
25 0. es OF O35 
36 0.30 0. O35 
49 0.38 OF 0.69 
64 0.47 0. 0.70 
8] O55 0: OS 
100 0.63 Ox 0.79 
12] Ol Al OQ. 0.83 
144 Or 7 Q. 0.93 
169 0.82 OQ. 1.00 
196 O ted Q. 1.00 
CE5 0.90 OQ. 0295 
250 O¢9 3 QO. 100 
289 0.96 Oz 0.99 
324 0.97 ve 1.00 
361 0.98 Sle e018 
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Figure 3.3: Probability of Detection vs. Number of Sonobuoys 
for a Figure of Merit of 65 
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Table 3.3: Probability of Detection vs. Number of Sonobuoys 
for a Figure of Merit of 65 


Number of Probability of Detection 
Sonobuoys 
Analytical Model B 

4 0.06 QO. 0.16 

9 a Ge: QO. Oe 218. 
16 Or! OF 0.36 
25 erst OQ. 0.54 
36 0.42 OQ. Oo 
49 Oss 0. 0.84 
64 O65 Ox 0G 
8] Oy OF ee 
100 0.80 an 1.00 
2) 0.86 On Ie (08 
144 0.9] 4). 1.00 
169 0.94 0. Go 
196 0.97 QO. 1.00 
Zap 0.98 0. 1.00 
256 0.99 0. 1.00 
289 1.00 QO. 1.00 
324 00 On 1.00 
S15) 1.00 OF 1.00 
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Figure 3.4: Probability of Detection vs. Number of Sonobuoys 
for a Figure of Merit of 70 
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Table 3.4: Probability of Detection vs. Number of Sonobuoys 
for a Figure of Merit of 70 


Number of Probability of Detection 
Sonobuoys 

Analytical Model B 

4 0.09 oF 0.16 

9 0.20 oF 0.4] 
16 Cue: 0. 0.70 
25 0.48 0. 0.95 
36 O76] OF 1.00 
49 0273 0. Cle 
64 0.84 Q. 60 
8] 0.90 OF 1.00 
100 0.95 QO. 1.00 
al 0.98 QO. 1-60 
144 0.99 oF 1.00 
169 1.00 0. 1.00 
196 1 #60 0. ato) 
aa 1 200 OF 1.00 
256 1900 0. 1.00 
289 1.00 0. ue 
324 1 200 [2 1.00 
361 1.00 QO. 1.00 
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Figure 3.5: Probability of Detection vs. Number of Sonobuoys 
Sensitivity Analysis for a Figure of Merit 60. L, S, and B, 
represent o,, 0., and o, respectively 
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Table 3.5: Probability of Detection vs. Number of Sonobuoys 
Sensitivity Analysis for a Figure of Merit of 60 


Number of Probability of Detection 
Sonobuoys 

G, =O. 0, = 2./ 0 =2./ of = 4.6 

Cee cet 6«6Oe 7 CC O.=2¢./ Gd, = 4.6 

O,=2./ og =2./ 0, =/7.0 og = 4.6 
4 0.09 0.08 0.08 0.08 
9 0.09 0.1] Or] 0.09 
16 OS 0.18 0.14 0.14 
2s, Og ote 0.30 0.26 O30 
36 0.29 0.34 O27 OR (28 
49 0.45 0.42 0.29 0.43 
64 0.45 0.5] 0.47 0.54 
8] Oe sre! 0.64 0.59 0.63 
100 0.56 0.69 0.68 O57 
121 0.64 On| 0.60 0.59 
144 0.65 O96 0.70 0.69 
169 OG 0.80 0.76 Oy 
196 0.69 0.86 0.83 0.84 
225 O60 076 0.79 02/9 
256 0.74 0.9] 0.91 0.88 
289 ays 0.87 0.84 0.83 
324 ONG 0.90 0.89 0.84 
36] Umes 0.93 0.95 0.9] 


re rt Erp || 
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IV. SUMMARY 


The purpose of this thesis is to investigate the effect of using the 
lambda-sigma jump model] in the detection component of APAIR. Two of the 
models that were used in the investigation are computer simulations. The 
simulation model referred to as Model A simulates the acoustic detection 
process similar to the way the Navy’s APAIR model does, whereas the 
Simulation model referred to as Model B differs from this in that, in 
addition, it incorporates a lambda-sigma jump process. The third model 
that was used is an analytical model called the random search model. This 
model provided a check on the structural accuracy of the two simulation 
models. 

As the FOM decreases the number of sonobuoys required to detect the 
submarine increases, for a given probability of detection. This was true 
for the two simulation models and the analytical model and for every value 
of FOM investigated, the relationships between the three models remained 
the same. Model B consistently produced higher probability of detection 
values for a given number of sonobuoys than did the other two models. 
Model A however, consistently produced lower probability of detection 
values and the analytical model produced values falling between the two 
Simulation models, for a given number of sonobuoys. 

The differences in the data generated by Model A and Model B were 
graphically significant. This indicates that the way Model A simulates 


Signa] excess fluctuation affects the model’s estimate of the number of 
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sonobuoys needed for a given probability of detection. The estimate of 
the number of sonobuoys needed for a given probability of detection using 
Model B, the lambda-sigma jump process model, was consistently less than 
the estimate obtained using Model A. 

In Model A, changing the values of the three standard deviation terms, 
0., 0, and og, significantly affected the models estimate of the number of 
sonobuoys required to gain a detection for a given probability. 

Finally in closing, for this scenario it seems the analytical model 
might be as capable and certainly less costly with respect to the two 
Simulation models in estimating the number of sonobuoys required to gain 
a detection. Model B consistently produced higher values of P, than 
Model A. This could be due to the fact that with Model B there is greater 
independence between successive values of X(t) whereas with Model A that 
is not the case. As the independence between successive values of X(t) 


increases so should the probability of detection. 
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APPENDIX A: 
THE SONAR EQUATION 


One of the basic forms for the passive sonar equation is: 

(A.1) SE = 10 log(S/N) - RD 
where: 

(A.2) 10 tog(S/N) = SL - TL - (AN - DI) 
and: 

SL = source level 

TL = transmission loss 

AN = ambient noise 

DI = directivity index 

RD = recognition differential 

SE = signal excess 
The unit of measure for the terms in the sonar equation is the decibel 
(dB). 

As can be seen from Equation A.1, signal excess depends on the 
following factors: source level; transmission loss; ambient noise; 
directivity index; and recognition differential. Source level is a 
measure of the amount of acoustic power that is radiated by a signal. 
Transmission loss is a measure of the energy loss that occurs when the 
sound travels from the source to the receiver. Ambient noise refers to 
all the acoustic energy that is received by a receiver that is in the 
ocean environment that is not part of the signal. Directivity index is 
a measure of how well the receiver can discriminate against noise arriving 


from directions other than that of the signal source. Recognition 
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differential is the signal-to-noise ratio in decibels at which the 
operator will detect the signal 50 percent of the time at some specified 
false alarm rate. In other words, it is the effectiveness of the sonar 
operator and the sonar equipment in detecting the signal. Signal excess 
is the difference between the signal-to-noise ratio in decibels and the 
recognition differential as seen in equation A.1. Frequency of merit, 
FOM, is defined as follows: 
(A.3) SE = FOM - TL. 

It is the transmission loss for which the signal-to-noise ratio in 
decibels is equal to the recognition differential. FOM is used with 


transmission loss curves in estimating sonar performance. 
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APPENDIX B: 
SENSITIVITY ANALYSIS 


To statistically determine if the values of o,, o, and o, have an 
effect on the probability of detection, an analysis of independence was 
done using a contingency table [Ref. 9:pp. 323-332]. The data used for 
this analysis, which is presented in Table B.1, was generated using Mode] 
A with the number of sonobuoys equal to 100 and an FOM equal to 70 dB. 
The contingency table is presented in Table B.2 and the X* statistic 
generated form this data is 118.6802 with 3 degrees of freedom. The 
results indicate that at « = 0.05 the null hypothesis, which is the 
probability of detection is independent of the values used for the 


Standard deviation terms, cannot be accepted. 
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Table B.l: Probability of Detection vs. Number of Sonobuoys 
Sensitivity Analysis with the Number of Sonobuoys = 100 


and FOM = 70 
Simulation Probability of Detection 
Run Number 
= 7.0 oO, = 2.) 0 =2.)/ oO, = 4.6 
=2.] oO, = 7.0 O0.= 2.) oO, = 4.6 
=2./ oO, =2./ 0o,=/7.0 oo, = 4.6 


oOooo0°oo 0 0 0 © 
qoooooooo eo © 


OoOwo onan on PP WP 
Oo ea 2 2 © 02 a@ @ © © 


— 
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Table B.2: Contingency Table of Probability of Detection with 
Number of Sonobuoys = 100 and FOM = 70 


ao=7.0 a =2.) oF =2.7) ao = 4.6 
Q=2.) oF = 7.0 op = 2.) GO, = 4.6 
On, = 2./ G6, = 2./ ao, = 7.0 On, = 4.6 

Py 838 909 894 874 

ee 162 091 106 126 


34 


APPENDIX C: 
PROGRAM LISTING 


Ke KK KK KK HK KK KK KK KK KKK KK KKK KKK KK KKK KKK KKK KKK KKK KK KKK KK KKK KKKKKEKKEKKEKER 


= rnis 
* sono 


program computes the probability of detection of a 
buoy field against a single submarine. 


* 
* 


KH He Ke He He He Ke Ke Ke HEH He He HK HK KKH KKH KKK KKH KKK KEKKKKKEKKKE KKK KEKKKEKREKKKKEKKKEKKKEKEE 


terete ett 


+ 
+ 
+ 


C INI 


PROGRAM DETECT 


REAL NTLV,SUBSPD, VSPD, XBNDRY, YBNDRY, XDAT, YDAT, 

STDV1, STDV2, STDV3 , SUBH, VARSH, SCALE, X1(1), 
SSL(1,15),AN(1),RD(1),DI(1),X2(1),X3(361), 
BLIFE(361),1TL(1,20), TBUOY(361) , SUM(361,1,361),PV3(3,3), 
XBUOY (361), YBUOY(361),DEC(6), XSUB, YSUB, TLINT, 

SHEAD, TIME, F, DIST, SBSPD,RSUB, SSD, SLD, 

SD, TLOSS, SPL, SPD,SQ,C,DSUB, DEGRAD, SL(361), 
FS1,RI,SP1,FOM,ISL,PDET, SAMPLE, 

RAND] (1) ,RAND3(1),RAND4(1),RANDS(1) 


PEGE ewes la MN LAMB, INTVL,1S1,1S2,1PRINT,NIVET, IMAX, 
Peo Omen tot, )S01,1S02,1SP1,1SP2,1BY, IDUMSNT, 
MAXINT, IDET, IZERO,NLINE,NFREQ,NREPS,NBUOY, ITER, ISTEP, 

JE BRE OE 


TIALIZE VARIABLES. 


DATA DSEED/123456713/ 
IDET =0 

DEGRAD =3.14159/180. 
J =] 

IZERO =0 


C READ IN DATA. 


+ 


C CAL 


READ(50,*) NLINE,NFREQ,NREPS, ITIME,IDEL,SAMPLE,NTLV, TLINT 
READ(50,*) SUBSPD,SCALE, XBNDRY, YBNDRY, ISEED, XDAT, YDAT 
READ(50,*) STDV1,STDX1,STDV3,MAXINT, IPRINT,VSPD 
READ(50,*) SUBH, VARSH, IMAX 

READ(50,*) ((SSL(I,J), J=1,15), 1=1,NFREQ) 

READ(50,*) (AN(J), J=1,NFREQ), (RD(J),J=1,NFREQ), 

(DI(J),J=1, NFREQ) 

READ(50,*) (DEC(J),J= 

READ(50,*) ((TL(I,d "5 ade vert) 


CULATE AND OUTPUT FOM. 
FOM = SSL(1,1) - AN(1) + DI(1) - RD(1) 
WRITE (45,290) FOM 
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290 FORMAT (2X,’FOM = ’,F6.2) 
WRITE(45, 300) 
300 FORMAT (2X,’NUMBER OF BUOYS’ ,4X,’PROBABILITY OF DETECTION’) 


C LAY IN SONOBUOY FIELD 


DO 10 L = 2, IMAX 
NBUOY = L**2 
IDET = 0 
DO 20 I = 1,NBUOY 
DO 30 J = 1,NFREQ 
ILINE(I,J) = 0.0 
30 CONTINUE 
20 CONTINUE 
SCALE = 2*YBNDRY/L 
DO 40 M = 1,L 
DO 50 N= 1,L 
1 = ((M-1)*L+4N) 
XBUOY ( (M-1)*L+N)=(XDAT-XBNDRY-SCALE/2)+SCALE*M 
YBUOY ( (M- 1) +L+N) = (YDAT- YBNDRY-SCALE/2)+SCALE*N 
TBUOY((M-1)*L+N)= 0.0 
BLIFE((M-1)*L+N)= 1500.0 
50 CONTINUE 
40 CONTINUE 


DO 60 I =1,NBUOY 
TBUOY(I) = 0.0 
60 CONTINUE 


CINE AL TZ Ess 
ISEED=ISEED + DSEED 
C PERFORM SIMULATION FOR THE NUMBER OF REPLICATIONS INDICATED 


DO 70 ITER=1,NREPS 
WRITE(6,310) ITER 
310 FORMAT(’ ’,’ ITERATION NUM=’ , 13) 


C INITIALIZATION OF VARIABLES 


INTVL=INT (SAMPLE) 
ISl=1 
IS2=IS1+INTVL 
CALL LRND(7,RAND1(1),1,1) 
SBSPD=SUBSPD+2*(RAND1(1)-0.5)*VSPD 
SLD=-1.0 
DO 80 I =1,NBUOY 
NIVET=0 
DO 80 J =1,NFREQ 
ILINE(I,J)=0 
iFLAG=0 
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80 

320 

900 

C FIND 


100 
C FIND 


110 


PV3(J,1)=ERG(STDX1, ISEED) 
PV3(J,2)=ERG(STDX1, ISEED) 
DO 80 K “I, MAXINT 
SUM(I, J,K)=-250 
CONTINUE 
IF(IPRINT.NE.1) GO TO 500 
WRITE(40, 320) 
FORMAT(’ ’,//,3X, ITER’ ,3X, TIME’, 2X, ‘BUOY’ ) 
CONTINUE 


LONG TERM GEOGRAPHIC FLUCTUATIONS. 

DO 100 J=1,NFREQ 
X1(J)=ERG(STDV1, ISEED) 

CONTINUE 

SONOBUOY PERFORMANCE. 

DO 110 I=1,NBUOY 


X3(1)=ERG(STDV3, ISEED) 
CONTINUE 


C GENERATE SUBMARINE LOCATION AND HEADING. 


330 


CALL LRND(9,RAND3(1),1,1) 
CALL LRND(3,RAND4(1),1,1) 
CALL LRND(4,RAND5(1),1,1) 
XSUB=XDAT+2 .0*(RAND3(1)-0.5 

(1)-0. 


l, 
l, 
1, 
Ws 
1) 
(1 


5) *XBNDRY 
YSUB=YDAT+2.0*(RAND4(1)-0.5)*YBNDRY 
SHEAD=(2*VARSH*(RAND5(1)-0.5)+SUBH) *DEGRAD 


IF (IPRINT.NE.1) GO TO 520 
WRITE(40,330) XSUB, YSUB, SHEAD, SBSPD 
FORMAT(4(4X,F6.2)) 


FOR EACH REPLICATION IN THE SIMULATION, TIME STEP THROUGH THE 


TIMES WHEN SIGNAL EXCESS IS COMPUTED FOR EACH SONOBUOY, TO 


C 
C ‘ENCOUNTER’. THE TIME STEPS, OR THE SAMPLING INTERVALS, ARE 
C 
5 


DETERMINE IF A DETECTION HAS OCCURRED. 


920 


DO 120 ISTEP=1,ITIME,IDEL 


C FIND SHORT TERM ENGAGEMENT FLUCTUATIONS. 


540 


TIME=FLOAT( ISTEP) 

DO 130 J=1,NFREQ 
IF(ISTEP.LT.1IS2) GO TO 530 
IF(J.LT.NFREQ) GO TO 540 
1S]=1S2 
IS2=IS2+INTVL 
PV3(J,1)=PV3(J,2) 
PV3(J,2)=ERG(STDX1, ISEED) 


Su) 


530 FS1=1S1 
IF(SAMPLE.£Q.0) WRITE(6,*) ‘ERROR SAMPLE’ 
F=(TIME-FS1)/SAMPLE 
X2(J)=(1.0-F}*PV3(J, 1) 4F*PV3(d. 2) 

130 CONTINUE 


C DETERMINE NEW SUBMARINE LOCATION AFTER TIME STEP INCREMENT. 


DIST=(SBSPD/60. )*IDEL 
X SUB=X SUB+ (DIST*SIN(SHEAD) ) 
YSUB=YSUB+(DIST*COS (SHEAD) ) 


C COMPUTE THE TRANSMISSION LOSS FOR EACH SONOBUOY. 


DO 140 I=1,NBUOY 
IF(TIME.LT.TBUOY(I)) GO TO 140 
IF(TIME.GT.TBUOY(I)+BLIFE(I)) GO TO 140 
RSUB=SOQRT( (XSUB-XBUOY (I) )**2+( YSUB- YBUOY (1) )**2) 
IF(TLINT.£0.0) WRITE(6,*) ’ERROR TLINT’ 
SSD=RSUB/TLINT 
ISD1=INT(SSD) 


C DETERMINE IF THE SUBMARINE IS OUT OF RANGE OF THE TRANSMISSION 
C 6LLOSS: VAWUESe 


IF(ISD1.GE.NTLV) GO TO 140 
SD=SSD-FLOAT(ISD1) 
ISD2=ISD1+1 

IF(ISD1.£Q.0) ISD1=1 
NIVET=0 


C PERFORM COMPUTATIONS FOR EACH FREQUENCY. 


DO 150 J=1,NFREQ 
TLOSS=(TL(J, [SD2)-TL(J, [SDde )*SR+ TLE ( Je Sia} 


C COMPUTE SOURCE LEVEL. 
SP1=SBSPD/2.0 
ISP1=INT(SP1) 
ISP2=ISP1+] 
SPD=SP1-FLOAT(ISP1) 
SL(J)=(SSL(J, ISP2)-SSL(J, ISP1))*SPD+SSL(J, ISP1) 
SLD=SBSPD 

C COMPUTE SIGNAL EXCESS. 


SQ=SL(J)-AN(J)-TLOSS-RD(J)+DI(J)+X1 (J) +X2 (J) +X3 (1) 
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C INTEGRATE SIGNAL EXCESS. 


IF(SQ.LT.0.0) GO TO 550 
IF(ILINE(I,J).EQ.J) GO TO 150 
IBY=ISTEP-INT(TBUOY(I))+1 
IDUM=MINO( IBY, MAXINT) 
IF(IDUM.LT.2) GO TO 560 
DO 160 K=2,IDUM 

C=FLOAT(K) 
IF(C.EQ.0) WRITE(6,*) ‘ERROR C’ 
SUM(I,J,K)=((C-1.0)*SUM(I,J,K-1)+SQ)/C 
160 CONTINUE 
560 SUM(I,J,1)=SQ 
IFLAG=ILINE(I,J) 


C DETERMINE IF SUBMARINE IS DETECTED. 


DO 170 K=1,1DUM 
IF(SUM(I,J,K).GE.DEC(K)) ILINE(I,J)=d 
FOM = SQ + TLOSS 
170 CONTINUE 
eee) «Ede. ) Goer 5/0 
GO TO 140 
570 CONTINUE 
NI =+] 
IF(IFLAG.EQ.0.AND.IPRINT.EQ.1) 
+ WRITE(40,340) ITER, ISTEP,NI 
IF (IPRINT.EQ.1) WRITE (40,350) FOM 
350 FORMAT(2X, “FIGURE OF MERIT = ’,F25.9) 
IFLAG=ILINE(I,d) 
DO 180 K=1,NFREQ 
NIVET=NIVET+ILINE(I,K)/K 
180 CONTINUE 
IF(NIVET.EQ.NLINE) GO TO 580 
GO TO 140 


C IF THE FREQUENCY WHICH WAS LOST WAS HELD BY THAT BUOY CALL 
C “LOST CONTACT’ ON THAT BUOY. 


550 IF(ILINE(I,J).NE.J) GO TO 590 


Te@eeRINT.&0.1) WRITE(40,340) ITER, ISTEP,NI 


C DECREASE THE NUMBER OF LINES HOLDING CONTACT BY THAT BUOY BY 
C ONE. 


NIVET=NIVET-1 


IFLAG=0 

590 CONTINUE 
ILINE(I,J)=0 

150 CONTINUE 
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140 CONTINUE 
C CONVERT SUBMARINE HEADING INTO DEGREES. 


DSUB=SHEAD*(180./3.14159) 
IF(IPRINT.NE.1) GO TO 120 
IF(ISTEP.GE.600) GO TO 120 
IF (ISTEP.EQ.1S1)WRITE(40, 340) ITER, ISTEP,NI, 
+  RSUB,XSUB, YSUB, DSUB 
120 CONTINUE 
GO TO 600 


C COMPUTE RANGE OF THE SUBMARINE TO THE SONOBUOY. 


980 RSUB=SQRT ( (XSUB-XBUOY (1) )**2+(YSUB-YBUOY (1) )**2) 
TIME=FLOAT (ISTEP) 


C OUTPUT SUBMARINE SPEED, POSITION, AND HEADING. 


IF(IPRINT.EQ.1) WRITE(40,360) SBSPD,XSUB,YSUB,DSUB 
360 FORMAT ("’, "SUBMARINE SPEED=",F4.1,2X, ‘XSUB=’, 
+ F9.2,2x,’YSUB=’ ,F6.2,2X, ‘DSUB=" ,F6.2) 


C OUTPUT RANGE OF SUBMARINE TO SONOBUOY. 


IF (IPRINT.£Q.1) WRITE (40,370) RSUB 
370 FORMAT(’ ’,’DISTANCE BUOY TO SUB=’ ,F5.2) 


C OUTPUT TIME OF DETECTION. 


IF (IPRINT.£Q.1) WRITE (60,380) ISTEP 
380 FORMAT(’ ’,’DETECTION AT TIME=’, 13) 
IDET=IDET+1 
600 CONTINUE 
70 CONTINUE 
IF(NREPS.LE.0) GO TO 610 


C COMPUTE AND OUTPUT PROBABILITY OF DETECTION FOR THE 
C CORRESPONDING SONOBUOY FIELD SIZE. 


PDET=FLOAT( IDET) /FLOAT(NREPS) 
610 CONTINUE 
WRITE(45,390) NBUOY,PDET 
390  FORMAT(9X,13,20X,F6.3) 
340  FORMAT(’ ’,14,4X,14,3X,13,4X,F6.2,6X,F6.2,5X,F6.2,4X, F6.2) 
10 CONTINUE 
STOP 
END 
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C THIS FUNCTION COMPUTES A STANDARD NORMAL RANDOM VARIATE. 


700 


FUNCTION ERG(SIG, ISEED) 
INTEGER LI 
REAL SUM, SIG, RAND6(1) 
SUM=0 
DO 700 Ll=1,6 
CALL LRND (ISEED,RAND6(1),1,1) 
SUM=SUM+RAND6 (1) 
CONTINUE 
ERG=1.41421*SIG*(SUM-3) 
RETURN 
END 


4] 


APPENDIX D: 
KEY VARIABLES 


Key variables of the FORTRAN program. 


AN(J) 
BLIFE(i) 
DEC(j) 
DI(j) 
IDEL 
IPRINT 
ISEED 
ITIME 
MAX INT 
NBUOY 
NFREQ 
NLINE 
NREPS 
NTLV 
RD(j) 
SAMPLE 


SCALE 
SSL(i, 3) 
STDV1 
STDV2 


ambient noise (dB) 

sonobuoy life (minutes) 

signal excess decrements (dB) 

directivity index (dB) 

Signal excess sampling interval (minutes) 
print output option 

random number initial seed 

Simulation time naees 

maximum number of samples per integration period (minutes) 
number of sonobuoys 

number of frequencies 

number of frequencies required for detection 
number of replications 

number of transmission loss values 
recognition differential (dB) 


interval for making random draws of short term fluctuation 
(minutes) 


scale factor for sonobuoy position 
submarine source level (dB) 
Standard deviation of long term fluctuation (dB) 


Standard deviation of short term fluctuation (dB) 
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STDV3 
SUBSPD 
SUBH 
TBUOY 
mei, 3) 
TLINT 
X1(j) 
X2(j) 
X3(i) 
XBNDRY 
XBUOY (1) 
XDAT 
YBNDRY 
YBUOY (1) 
YDAT 
VARSH 


Standard deviation of buoy to buoy fluctuation (dB) 


submarine speed (knots) 


Submarine heading (degrees) 


sonobuoy activation time (minutes) 


transmission loss (dB) 


transmission loss interval (nautical miles) 


long term fluctuation (dB) 


short term fluctuation (dB) 


buoy to buoy 


X 
X 


X 


y 
7 
yf 


coordinate 
coordinate 
coordinate 
coordinate 
coordinate 


coordinate 


fluctuation (dB) 
of spa boundary 
of sonobuoy 

of datum 

of spa boundary 
of sonobuoy 


of datum 


Standard deviation of submarine heading (degrees) 
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